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ABSTRACT 

We report on a re-analysis of archival data from the Very Large Array for a sample of ten long 
duration radio transients reported by Bower and others. These transients have an implied all-sky rate 
that would make them the most common radio transient in the sky and yet most have no quiescent 
counterparts at other wavelengths and therefore no known progenitor (other than Galactic neutron 
stars). We find that more than half of these transients are due to rare data artifacts. The remaining 
sources have lower signal-to-noise ratio (SNR) than initially reported by 1 to 1.5 a. This lowering 
of SNR matters greatly since the sources are at the threshold. We are unable to decisively account 
for the differences. By two orthogonal criteria one source appears to be a good detection. Thus the 
rate of long duration radio transients without optical counterparts is, at best, comparable to that of 
the class of recently discovered Swift J1644+57 nuclear radio transients. We revisit the known and 
expected classes of long duration radio transients and conclude that the dynamic radio sky remains 
a rich area for further exploration. Informed by the experience of past searches for radio transients, 
we suggest that future surveys pay closer attention to rare data errors and ensure that a wealth of 
sensitive multi-wavelength data be available in advance of the radio observations and that the radio 
searches should have assured follow-up resources. 
Subject headings: catalogs — radio continuum: galaxies — surveys 



1. INTRODUCTION 

A century ago, the study of variable stars was a lead- 
ing area of astronomy. With the increasing availability 
of large format optical detectors and inexpensive high 
speed computing this sub-field, as witnessed by projects 
such as AS AS, OGLE, Catalina Sky Survey, the Palo- 
mar Transient Factory, PanSTARSS and SkyMapper, is 
making a come-back. Radio astronomy appears to be 
poised for a similar growth. At meter wavelengths, com- 
mercially available signal processing chips make it fea- 
sible to image the entire primary beam of a dipole or 
a cluster of dipoles. The se technological innova tions lie 
at the heart of LOFA R (iRottgering et al.l 120031) MWA 
(|Lonsdale et alj 120091) and LWA (|Ellingson et aL|[2009h . 
At centimeter wavelengths the "large number, small di- 
ameter" (LNSD) array approach (made possible by in- 
expensive signal processing, advances in commercial RF 
technology, innovative ideas in the design of small diam- 
eter telescopes and phased array focal planes) has now 
been demonstrated to be a cost effecti ve method of build- 
ing high speed mapping m achi nes dWelch et alJ 120091 : 
iDewdnev et alJl200aiJona3l25o9LIOosterloo et alH l2009)l 
The LNSD approach has motiva ted a new generation o f 
radio facili ties: Apertif/WSR T (lOosterloo et alJUOOl): 
MeerKAT (|Booth et all 120091) : ASKAP (jJohnston et al.1 
120081) . 

We divide radio transients into four categories based 
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on two attributes. The first is the duration of the ba- 
sic phenomenon (shorter than or greater than a few sec- 
onds). The second is their location (within the Galaxy 
or extra-galactic). Roughly speaking the duration maps 
to coherent versus incoherent emission and the location 
to repeated versus cataclysmic events. 

Pulsars and related phenomenon (giant pulses, nulling 
pulsars, erratic pulsars, rotating radio transients, and 
magnetars) are the dominant category of short duration 
radio transients at meter and centimeter wavelengths. 
There are no secure examples of short duration radio 
transients that are located beyond the local Group. Flare 
stars and associated phenomena are prime examples of 
long duration radio transients of Galactic origin. 

The focus of this paper is long duration transients of 
extra-galact ic origin. Known e xamples in this group are 
supernovae (iWeiler et al lf2010h and gamma-ray burst af- 
terglows (jGehrels et al.l 120091 ) . In both cases, the radio 
emission arises as the fast moving debris interacts with 
the circumstellar matter. In Table [T] we summarize the 
areal density of radio-emitting supernovae (including the 
sub-classes) and GRB afterglows. Note the areal density 
of "live transients" (transients present at any given in- 
stant of time) of both supernovae and GRB afterglows is 
less than 0.05 per square degree. 

In 2007, Bower et al. [hereafter, B07] reported on the 
analysis of a single field observed every week as a part of 
the Very Large Array (VLA) calibration program. The 
observations were conducted at 4.8 GHz and 8.4 GHz and 
lasted 22 years. The 944 epochs and the weekly cadence 
makes this data set a most valuable set to probe the 
decimeter band for long duration transients at the sub- 
milliJansky level. These authors reported the discovery 
of eight transients found in only one epoch (hereafter 
"single-epoch"; duration, 20 minutes < t^ur < lweek) 
and two transients found in rolling 2-month searches 
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TABLE 1 

Long Duration Transient Populations 



Class 


Rise 


Decay 


D 


Host Rate 


Ref. 






yr 




mag deg -2 





Type II SNe 


0.1-1 


10 


100 Mpc 


16 


0.04 




Type Ib/c SNe 


0.1 


0.3 


50Mpc 


14.5 


5 x 10- 


-6 


SN1998bw-like 


0.1 


0.1 


300 Mpc 


18.4 


3 x 10~ 


-4 


SwJ1644+57-like 


0.1 


1 


z ~ 1.8 


21.7 


0.1 




Orphan Afterglows 


1 


1 


lGpc 


21.0 


10~ 2 




NS-NS mergers 


0.1-1 


0.1-3 


800 Mpc 


20.5 


5 x 10- 


-3 



Gal- Yam et al. (2006) 
Bereer et al. (2003) 
Soderberg et al (20061 
Zauderer et al. (2011) 
LevinsonetaL 
Nakar fc Piran 



(2002) 
' (2011) 



Note. — Dctcctability distance and rates have been calculated assuming a single snapshot at a 
flux density threshold of 0.3 mjy. See section[5]for details. D is the distance at which the typical 
transient will have a specific flux of 0.3 mjy. Host is the apparent magnitude of a galaxy with —19 
absolute magnitude at distance D. 



(hereafter, "multi-epoch" sources). 

Deep observations towards these sources were under- 
taken at optical, near-IR and X-ray bands. The most re- 
markable feature of the B07 sources is an absence of opti- 
cal and near-IR c ounterparts, desp ite deep sea r ches (B 07; 
IQfek et al.l l20ToT) . As noted bv IQfek et afl (l2010h all 
extra-galactic transients (regardless of the band at which 
the transient was discovered) have detectable optical 
counterparts, namely, their host galaxies. Furthermore, 
remarkably the areal density of transients live at any 
given time was estimated to be 1.5deg -2 (S > 0.37 mjy). 
This density exceeds that of all other known radio tran- 
sient source popula tions by an order of magnitude (or 
more); see Table [T] IQfek et al.l (|2010f) thus argued that 
the absence of an optical counterpart means that B07 
transients have to be repeating sources of Galactic ori- 
gin, and proposed that B07 transients are old neutron 
stars (which naturally satisfy the requirement of being 
optically almost invisible). 

Given that the search for transient and strong vari- 
ables is one of the primary motivators for the next gen- 
eration radio facilities (described earlier) it is important 
to critically investigate the B07 transients since this class 
nominally dominates over all other known classes of ra- 
dio transients (see Tabled]). To this end, here we report 
on a re-analysis of the original d ata of B07 (jjH Q and 
revisit the transient reported bv lOfek et alJl201lT T ^3.8|) . 
In $5] we present an update of the expected rate of ra- 
dio transients. In Sj6]we discuss these rates in relation to 
future synoptic radio imaging surveys and conclude. 

2. OBSERVATIONS AND RE-ANALYSIS 

The data used by B07 arose from a calibrator pro- 
gram carried out during the period 1983-2005. All ob- 
servations were made in the standard continuum mode 
with 100 MHz of total bandwidth in each of two adjacent 
50-MHz bands (IFs) at center frequencies of 5 GHz and 
8.4 GHz and in both hands of circular polarization. Sec 
B07 for more details about the full data-set. 

For the re-analysis we selected, from the archive, only 
the raw data relevant to the transients reported in B07. 
This means the eight epochs from which the single-epoch 
transients were first found and the 3+8 data sets from 
which the two multi-epoch transients were found. Data 
were taken at other radio frequencies in about half of the 
cases. Some details of the single epoch and two-month 
transients can be found in Table [2]) . 



For the re-analysis we used AIPS3 (|Greisenll2003l) . The 
data reduction and imaging followed the same path used 
by B07 with a small exception. B07 employed AIPS for 
the flagging and analysis of the single-epoch t ransients, 
and used the Miriad package (|Sault et ali i 995) for imag- 
ing the two-month averages. We endeavored to make 
the calibration and the flagging of UV data (AIPS task 
TVFLG) data for each epoch as uniform way. Following 
these steps we ran each raw visibility data set through 
the VLA pipeline (VLARUN). 

No flux density calibrator was observed during any 
epoch of these test observations. Following B07, the flux 
density of the phase calibrator (B1803+784) was fixed to 
be2.2Jy (5 GHz) and 2.8 Jy (8.4 GHz). For those epochs 
with 22 GHz and 1.4 GHz observations the flux density 
of the phase calibrator was taken to be 3 Jy and 2 Jy, re- 
spectively. It is evident from the strong variations in the 
radio light curves for B1803+7843 that these mean val- 
ues are only approximate. Our reinvestigation confirm 
that at least during the period 1981-1999 the variation 
was less than 15%. Fortunately, an accurate flux density 
scale is not crucial for our analysis since we report results 
in terms of the signal-to-noise. 

Following B07 we applied a Gaussian weighting to the 
visibility data in order to limit the effects of bandwidth 
smearing. This was done by applying a 150-kA taper 
to all visibility data prior to imaging (IMAGR). For each 
field we required that a source be present in both fre- 
quency bands (IFs) with similar flux densities and with 
similar positions. Images were made with extra large 
fields-of-view. The wide field-of-view is necessary to re- 
duce the effect of side-lobes that can mimic sources in 
narrow fields. These final analysis images had a size 
of about 40-arcmin at 5 GHz data and 27-arcmin at 
8.4 GHz. For guidance, the full-width at half power for 
VLA antennas is 45-arcmin/j/(GHz), or 9.3-arcmin at 
5 GHz and 4.3-arcmin for 8.4 GHz. Measurements of the 
V LA beam power response be yond the first null are given 
in I Cotton and Per lev I (|2010f ). while the polynomial co- 
efficients needed to correct for the primary beam attenu- 
ation can be found in the AIPS task PBCDR. We do note 
that these corrections are uncertain for large angular dis- 
tances from the beam center. All data were taken in the 
B1950 coordinate system. We stayed in the B1950 sys- 

6 http://www.aips.nrao.edu/ 

7 U. Michigan Radio Astronomy Observatory database 
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tern throughout calibration and imaging. 

3. FINDINGS 

Below we offer a detailed report for each of the ten 
sources reported in B07. Summarizing our results for the 
impatient reader we find four of the eight single-epoch 
transients [RT 19840502 f ^TT]) , RT 19860122 (g33J), 
RT 19840613 (gSHJ, RT 19860115 (gp])] and one of the 
two-month transients (RT 19870422; ^3~5l) to be artifacts. 
Our re-analysis finds that the remaining long-duration 
transient, RT 20010331 ( £)3.6I) did not have a significant 
detection when imaged with AIPS (SNR of 3.2) but the 
Miriad imaging yields an SNR of 7.4. 

For the remaining single epoch sources [RT 19920826, 
RT 19970205, RT 19970528, RT 19990504] our analysis 
(undertaken by DAF; see ^5] for a summary of our data 
reduction) finds reduced SNR. The discrepancy with re- 
spect to the results reported in B07 concerned us and so 
another author (GB) reanalyzed these fields. Some of the 
discrepancies arise from UV data flagging but even when 
the same flagging is used different SNR algorithms yield 
measures which differ by 1 to 1.5 a. Since the sources 
lie close to threshold of detection even a small shift of 
a single a has an exponential effect in their confidence. 
In effect, the reality of this group (taken as a whole) ap- 
parently depends on details of algorithms in AIPS and 
Miriad - the investigation of which is beyond the scope 
of this paper. Below in ^3.7.11 17.7.31 we detail the analy- 
ses and a summary of the SNRs of these sources can be 
found in Tabled 

We also re-investigated th e lone transient candidate 
reported in lOfek et al.l (|2011l) . Our re-examination which 
now correctly includes the number of independent beams 
that were searched for shows that there is a four percent 
probability that this candidate is due to noise (see £)3.8[) . 
In the next section (Q we synthesize these findings and 
present our conclusions about the B07 transients. 

3.1. RT19840502 

B07 report finding a transient close to the pointing 
center (13") with a primary beam-corrected flux density 
of 448 ± 74 /iJy, or a SNR=6.1. We imaged the calibrated 
data set and confirmed the presence of emission at this 
level at the reported position. However, we identified 
a phase center artifact in the visibility data for the left 
hand polarization of the upper IF band. The effect of this 
artifact in the image plane is to create strong positive and 
negative slidelobes, with the positive feature identified as 
a transient source (Figure [1]) . When the upper IF band 
data are removed the resulting image is noise-like and the 
peak flux density at the nominal position of the reported 
transient is 191 ± 97/iJy. Additional observations were 
taken during the same epoch at 15 GHz. The peak flux 
density at the same position is —243 ± 200 /iJy. 

3.2. RT19840613 

B07 report finding a transient coincident with a pos- 
sible host galaxy (z = 0.040) and with a primary beam- 
corrected flux density of 566 ± 81 fiJy, or a SNR=7.0. 
Our deconvolved image shows a source at that location. 
Gaussian fitting suggests that the source is resolved and 
this conclusion is supported with an integrated flux den- 
sity (715 ± 218 ^tJy) being clearly larger than the esti- 
mate of the peak flux density of 388 ± 82 fiJy. Imaging 
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Fig. 1.— Image of RT 19840502 marked by cross. The fea- 
ture at the center (0,0) is an artifact ("phase center"). Once 
the bad visibility data is removed the transient candidate is 
not visible on the final image. Contours are displayed in steps 
from — 1, —0.75, —0.5, 0.5, 1.0, 1.5 mjy, with negative (posi- 
tive) contours given by dashed (sold) lines. The synthesized 
beam size is shown in the lower left. 

the lower and upper IFs separately we find another dis- 
crepancy. The peak flux density in the lower IF is four 
times weaker than the upper IF band. 

The likely source of the problem can be seen by com- 
paring the dirty image with the dirty beam (Figure [5]) . 
RT 19840613 appears to be an uncleaned side-lobe of 
J150123+781806, one of the brighter persistent sources 
detected in 452 images made by B07 at 5 GHz. The puta- 
tive transient is 56" away from J150123+781806 to the 
northwest, close to a local maximum (10% of peak) in 
the dirty beam at this location. This side-lobe artifact is 
stronger in the upper band but it is still present in the 
lower band. 

Deconvolution does not fully remove the side-lobe from 
the image and the effect is to produce a false transient. 
We investigated whether the artifact is due to short- 
timescale (~ lOmin) variability of J150123+781806 but 
after dividing the visibility data in half (by time) and 
re-imaging, we found no evidence for variability. The ar- 
tifact may be due to some low level interference picked 
up by some antennas baselines but we were not able to 
identify the bad data. It is possible that RT 19840613 is 
a real transient that unfortunately lies on the side-lobe 
of the dirty beam, but its Gaussian fit parameters and 
its variation in the lower and upper sidebands suggest 
that it is not a real source. 

Additional observations were made during this epoch 
at 1.5 GHz. The data quality is good and the images have 
no obvious artifacts. No source is visible at the transient 
position. The peak flux density is 133 ± 113 /iJy. 

3.3. RT19860115 

B07 report a transient with a primary beam-corrected 
flux density of 370 ± 67/iJy, or a SNR=5.5. Some of 
the same issues with the image of the previously dis- 
cussed RT 19840613 were also seen for RT 19860115. In 
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TABLE 2 

Single Epoch and Two-Month Transient Candidates 



Name 


Type 


Freq. 


FWHM/2 


Ad 


Beam 




Candidate 




(GHz) 


(arcmin) 


(arcmin) 


(arcsec) 


Notes 


RT 19840502 


SE 


4.9 


4.6 


0.22 


6.0 


Phase center artifact. 


RT 19840613 


SE 


4.9 


4.6 


3.3 


5.7 


Side-lobe of bright source. 


RT 19860115 


SE 


4.9 


1.6 


1.3 


14.8 


Side-lobe of bright source. 


RT 19860122 


SE 


4.9 


4.6 


1.6 


14.5 


Artifact. Lower IF is bad. 


RT 19870422 


2M 


4.9 


4.6 


8.0 


12.8 


Artifact. Bad pointing. 


RT 20010331 


2M 


8.5 


2.6 


1.1 


1.5 


No detection (see i|3.6>. 


RT 19920826 


SE 


4.9 


1.6 


2.0 


21.2 


D SNR: 6.4; 5.8. 


RT 19970205 


SE 


8.5 


2.6 


4.1 


1.4 


B SNR: 7.7; 5.7 (CA) 


RT 19970528 


SE 


4.9 


4.6 


6.8 


3.9 


CnB SNR: 7.5; 5.6. (CA) 


RT 19990504 


SE 


4.9 


4.6 


8.9 


18.9 


D SNR: 7.3; 5.7. 


RT 19970528 


SE 


8.5 


2.6 


6.8 


1.3 


CnB No detection (BA) 


RT 19990504 


SE 


8.5 


2.6 


8.9 


8.3 


D No detection (BA) 


Note. — Starting 


from the 


left the columns are as 


follows. 


The name of the transient as 



RT YYYYMMDD where YYYY is the UT year, MM is the month index and DD is the day num- 
ber at which the transient was first detected; the type of transient: single epoch (SE) or two month 
average (2M); the center frequency in GHz; one half of the full-width at half-maximum of the primary 
response beam in arc minutes; the offset of the transient from the phase center, also in arc minutes; the 
beam size in arc seconds computed as geometric mean of the major and minor axes;. The last column 
report the array configuration, two SNRs (for sources detected and reported as such in B07) and some 
comments. The left SNR (in italics) are the SNR from B07 and the right SNR (in typewriter font) 
arc SNR resulting from the work presented here (sec [j3]for details). The comments are as follows: CA 
for severe Chromatic Aberration and BA for severe Beam Attenuation. The first group of transients 
are either artifacts or clearly very low SNR. Following that (separated by a line) are four transients 
which arc threshold sources. The last group is composed of higher frequency observations of two of 
these threshold sources and given the strong beam attenuation there was no expectation of detection. 
The entries were made merely for completeness. 



Figure [3] we show the dirty image along with the syn- 
thesized beam. The RT 19860115 appears to be an un- 
cleaned side-lobe of J150123+781806 and hence not a real 
transient. RT 19860115 lies at the same angular distance 
(3.3') and position angle (100° CCW CC) of a local max- 
imum in the side- lobe structure (25% of the peak beam). 
We were unable to identify the source of these strong 
side-lobes. As in the case of RT 19840613, we were able 
to rule out that the strong side-lobes originated from 
short-term variability of J150123+781806. 

3.4. RT19860122 

B07 report a transient with a primary beam-corrected 
flux density of 1586 ± 248 /zJy, or a SNR=6.4. Images 
made with all the visibility data do show emission at 
the position of RT 19860122 with the reported peak flux 
density. However, there does appear to be some erro- 
neous visibility data resulting in low-level rings in the 
image centered at the phase center. The bad data were 
traced to the lower IF. If the data for the lower IF band 
are removed there is no source at the reported transient 
position (peak flux density of 303 ± 189 /iJy). 

3.5. RT 19870422 

This is one of two transients found by binning in- 
dividual epochs into two-month averages. B07 report 
a 5 GHz transient with a primary beam-corrected flux 
density of 505 ± 83/iJy, or a SNR of 6.1. The source 
J150050+780945.5 is positionally coincident with a blue 
host galaxy (z = 0. 249) and has been suggested by 
iNakar fc Piranl (|2TJll as a candidate afterglow from a 
binary coalescence event. Three single epochs were used 
to form the average. They are in YYYYMMDD format: 
19870414, 19870422 and 19870429. A close inspection of 



the data reveals a problem. The last two of the three 
epochs used to find RT 19870422 were pointed to a dif- 
ferent part of the sky, an error that was made when the 
original observing files were written. The source is ac- 
tually J210133. 67+373528, and appears only as a tran- 
sient when added inappropriately to the data from the 
standard on-the-sky test field. In hindsight we recognize 
these to be the same type of errors as seen earlier (see 
lOfek et al.l[2TjT0l) . 

3.6. RT20010331 

This is the second of the transients which were found by 
binning individual epochs into two-month averages. B07 
report a 5 GHz transient with a primary beam-corrected 
flux density of 697 ± 94 //Jy, or a SNR=7.4. Separately, 
iCroft et all (f201lT) reported a marginal X-ray source at 
this position. 

Eight single epochs were used to form the average. 
They are in YYYYMMDD format: 20010306, 20010314, 
20010321, 20010328, 20010403, 20010411, 20010418 and 
20010425. We calibrated and imaged the eight epochs of 
observations from 2001 March 6 to 2001 April 25 (all in 
B-configuration) . While the four quiescent sources from 
B07 can be seen in this deep image, we do not iden- 
tify a significant source at the position of RT 20010331. 
The brightest peak within the synthesized beam (natu- 
ral weighting) is 42 ± 13 /iJy (or 3.2er). An independent 
reduction using the AIPS package by GB confirms the 
absence of this source. 

The first problem lies with the primary beam correc- 
tions reported in B07. RT 20010331 lies 4.4' from the 
phase center, close to the 10% response radius of the pri- 
mary beam. Our estimate for the rms noise of 13 /iJy 
lies within a few percent of the theoretical value. The 
beam-corrected rms noise in this case (at 10% response 
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Fig. 2.— Top: The dirty beam for the RT 19840613 field. 
This sub- image shows the northwestern side of the beam. The 
peak is located in the lower right hand corner. The contours 
are 5, 10, 15, 20, 25 and 30% of the peak flux. Bottom: 
The dirty image toward RT 19840613. A cross marks the 
location of transient. Contours are 0.13, 0.23, 0.33, 0.43, 
0.53, 0.63 mjy beam" 1 . Note that the transient candidate 
RT 19840613 lies at the same angle and position as a side- 
lobe from the bright source J150123+781806 located in the 
bottom right corner of this image. 



of the primary beam) would then be about 129 yuJy, not 
the 94/uJy given in B07. Rather than correcting at the 
10% radius, it appears that the flux density and noise in 
B07 were mistakenly corrected at the 14% power level. 
This multiplicative error has no impact on the signal-to- 
noise. 

In order to investigate this signal-to-noise discrepancy 
between our image and B07, we split the data in vari- 
ous ways (separate epochs, months of March and April 
epochs, adjacent IF bands) and re-imaged, looking for a 
bright peak. None were found. 

The discrepancy between B07 and the work reported 
apparently can be traced back to differences in the Miriad 
and AIPS imaging packages. Our calibrated visibility 
data, when processed through Miriad using nearly iden- 
tical imaging parameters as those in AIPS, gives a flux 
density of 91 ± 13 ^Jy (7a). We have no explanation for 
the discrepancy between the two results obtained from 
AIPS and Miriad. It is worth noting that the peak flux 
densities of the persistent sources identified by B07 agree 
in these images to 0.5a. For this paper we accept the 
analysis given here. 

There was one epoch (2001 March 6) in which data 
were also taken at 5 GHz. The peak flux density at the 
position of RT 20010331 is -27± 39/zJy. In summary, 
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Fig. 3.— Top: The dirty beam for the RT 19860115 field. 
This sub-image shows the southwestern side of the beam. 
The peak of the beam is located in the upper right hand 
corner. The contours are 7.5, 12.5, 17.5, and further in 
steps of 5% up to 37.5% of the peak flux. Bottom: Dirty 
image toward RT 19860115. A cross marks the location of 
the transient. The contours are 0.15, 0.25, 0.35, 0.45, 0.55, 
0.65 and 0.75 mjy beam -1 . Note that the transient candidate 
RT 19860115 lies at the same angle and position as a side-lobe 
from the bright source J150123+781806 located in the upper 
right corner of this image. 

we find no evidence to support that RT 20010331 is a 
significant detection. 

3.7. Remaining Four Single Epoch Sources 

3.7.1. RT19920826 

B07 report a transient with a primary beam-corrected 
flux density of 642 ± 101 /iJy, or a SNR=6.4. We con- 
firm a source at this position but with a slightly reduced 
SNR. Using natural weighting of the gridded visibilities, 
our measured peak flux density is 460 ± 80 /zJy, or 5.8o\ 
Despite the lower significance, there is some confidence 
that RT 19920826 is real since it appears in both IFs with 
comparable flux densities. Further investigation by GB 
shows that the difference in the SNR between B07 and 
the analysis here can be traced to differences in flagging 
of the UV data (of two specific antennas). 

3.7.2. RT 19970528 

B07 report a transient with a primary beam-corrected 
flux density of 1731 ± 232 ^Jy, or a SNR=7.5. These ob- 
servations were taken during a time when some antennas 
were being moved to the B configuration and so we ap- 
plied antenna position corrections (via AIPS task VLANT) 
before calibration. 

RT 19970528 is 6.8' from the center of the image, 
where the response of the antennas is only 16% of their 
peak. The uncorrected flux density is 270 ± 47. A point 
source search (AIPS task SAD) of the four million pix- 
els enclosed interior to a radius around this candidate, 
shows six other un-cataloged candidates with similar 
signal-to-noise. The primary-beam corrected flux density 
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1.7 ± 0.3 mJy is similar that from to B07. This estimate 
does not include a correction for temporal smearing due 
to the sources distance from the phase center, nor for 
the added uncertainty in the magnitude of the primary 
beam correction beyond the 20% point. Correcting for 
these terms we get 2.8 ± 0.5 mJy. 

GB re-investigated this source and found that with 
the B07 analysis the SNR varies with the approach 
(peak/rms, SAD, JMFIT) between 6.2, 7.4 and 7.1. (with- 
out any corrections discussed above). The corresponding 
SNRs for the image discussed above is 6.1, 5.1 and 6.8. 
The source is a bit "ratty" and this may explain the vari- 
ation in SNR. We therefore find this source to be a weak 
detection. 

Additional observations were made during this epoch 
at 8.5 GHz. The source lies close to the first null 
(6.4 ± 0.3' from the phase center) of the beam at this 
frequency. The attenuation by the primary beam is se- 
vere and hence the sensitivity is not sufficient to provide 
any strong spectral index constraints. 

3.7.3. RT 19990504 

B07 report a transient with a primary bea m-corrected 
flux d ensity of 7042 ±963 /xJy, or a SNR=7.3. ICroft et al.l 
(|2011f ) report an X-ray source in the vicinity of the radio 
source. Deep multi-wavelength data is consistent with 
the X-ray source arising from a QSO but located five 
arcseconds from the putative radio transient. 

We find the (uncorrected) peak flux density is 290 ± 
51 ^Jy (or an SNR of 5.7). RT 19990504 lies 8.9' from the 
phase center, close to the first null (11.1 ± 0.6' from the 
phase center) where we would not expect to find sources. 
The polynomial expressions used in AIPS to correct for 
the beam attenuation are increasingly inaccurate outside 
the 20% response radius, and they are not applicable 
close to the null. The azimuthally averaged measured 
value of the beam response is 1.8%, implying a flux den- 
sity correction factor of 55 x (|Cotton and Perlevl 120101) 
or 16.0 ± 2.8 mJy. 

As with RT 19970528, the analysis of B07 data by GB 
finds SNR of 5.6 (peak/rms), 7.3 (SAD) and 6.8 (JMFIT) 
and a similar variation with the image reported here. It 
is worth noting that the source may be extended and 
also that a visual inspection of the annular region re- 
veal a number of similar sources. We conclude that 
RT 19990504 is not a robust detection. 

Additional observations were made during this epoch 
at 8.5 GHz. An image was made from these higher fre- 
quency data but, as expected, no source was found. 

3.7.4. RT 19970205 

This is the only single-epoch transient identified 
at 8.5 GHz by B07. They report a primary beam- 
corrected flux density of 2234 ± 288 fiJy, or a SNR=7.8. 
RT 19970205 is 4.4' from the center of the image, where 
the response of the antennas is only 9.8% of their peak. 
There is an elongated source at this position with a peak 
flux density of 231 ± 41 /iJy (using natural weighting). 
There is some indication that the source is extended since 
the integrated flux density is about twice the peak flux 
density. 

The peak flux density was likely underestimated by 
B07 since they did not fully correct for bandwidth or 
temporal smearing effects. These data were taken in the 



BnA array with an integration time of 3 1/3 s and a 
50-MHz bandwidth. In this observing configuration, the 
effects of temporal smearing with this dump time will 
likely reduce the peak flux density by only a few percent. 
However, chromatic aberration is expected to be larger 
with the e ffect of smearing a point source along the radial 
direction (jPerlev et al.lll989D . 

Accounting for the bandwidth effect, we measure a 
peak flux density of 3.4 ± 0.6 mJy (SNR=5.7). Next, 
due to the uncertainty about the form of the primary 
response function beyond the 20% point (AIPS task PB- 
COR), we only know that the peak flux of RT 19970205 
lies in the range of 2.9-4.4 mJy. As noted above the 
source is likely extended and this may account for the 
variation in SNR (from 6 to 7.6). 

3.8. Other Surveys 

Ofek et akl f|201 lh presented a survey for 5 GHz radio 
transients at low Galactic latitudes. Most of the data 
were reduced in near-real time (2hr delay) and transient 
candidates were followed up with radio, visible light and 
X-ray instruments. The authors reported a single tran- 
sient candidate, J213622. 04 + 415920.3 with a peak spe- 
cific flux of 2.36 ± 0.41 mJy (5.8 a) and with no obvious 
optical counterpart. 

We have re-analyzed this data set and confirm that 
the SNR is 5.8. Considering that the search resulted 
from inspecting 1.1 x 10 7 independent beams ( CP|) and 
assuming Gaussian statistics, we find that the probability 
of the highest event found in these many independent 
beams is attributable to chance or noise is 3.6% (see fA"]l. 
Therefore, we advocate that this candidate is not a real 
event. 

4. THE REVISED B07 TRANSIENT RATE 

As noted in Sjl] the ten transients reported in B07 were 
remarkable for their apparent abundance (1.5 per square 
degree, at any given epoch) and the lack of quiescent op- 
tical counterparts. In $3] we found that five of the B07 
transients are artifacts arising out of (rare) data acqui- 
sition problems or imaging artifacts (see Table [2]) . We 
find that the 2-month transient RT 20010331 is not well 
detected. We arg ue that the single transient reported by 
I Ofek et ail (|201ll) is also not a robust detection. 

In $3] we present our analysis for the remaining 
four transients (all of which are single-epoch tran- 
sients): RT 19920826, RT 19970205, RT 19970508 and 
RT 19990504. We find SNRs which are typically 1-cr be- 
low that reported in B07. The revised and B07 SNRs are 
given in Table [2] We no te that an independe nt pipeline 
written in ParselTongue (jKettenis et al. 2006), a Python 
interface to AIPS, and run on three of these four source s 
confirms the lower SNR values found here (|Bell ll20TlT ): 
RT 19920826 (SNR=6.0), RT 19970528 (SNR=4.4), and 
RT 19990504 (SNR=4.0). 

The lowering of SNR (from between 7 and 8 to between 
5 and 6) has a pernicious effect when the number of inde- 
pendent beams which were search is included. In CdI we 
estimate this number to be n w 9 x 10 7 . In [JX]we derive 
the probability density function for the highest m values 
of n Gaussian random numbers^. In Figure|4]we plot the 

8 Theory informs us that the statistics of beam values or equiv- 
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density function for the highest value (m = 1) and the 
fourth highest value (m = 4). As can be seen from this 
Figure |U if the SNRs reported here are accepted then the 
global case for the remaining B07 transients is entirely 
weakened. If, on the other hand, the SNRs reported in 
B07 are accepted then the four transients reported in 
B07 do argue for a new class of radio transients. 
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Fig. 4. — The probability density function of the highest 
value (m = 1) and the fourth highest value (m = 4) of a 
population of n = 9 x 10 7 Gaussianly distributed random 
numbers with zero mean and unit variance. The probabil- 
ity density function plotted here is discussed in ijA] and the 
justification for n can be found in SjD] The dots represent 
the SNRs of the following single-epoch transients discussed 
in Tabled (from left to right): RT 19905054, RT 19920826, 
RT 19970205 and RT 19970228. (Top): SNRs as reported in 
this paper. (Bottom): SNRs as reported in B07. Note that 
the vertical location of the dots is arbitrary. For details of 
these density distributions see f|7\] 

The above approach of using a fixed threshold for all 
epochs does not result in optimal detection. In partic- 
ular, the threshold for a low resolution survey is lower 
than that for a higher resolution survey (since the lat- 
ter has a correspondingly larger number of synthesized 
beams). B07 addressed this problem by requiring that 
the probability of a false detection (PFD) in an individ- 
ual epoch was constant and less than N where N is the 
total number of images. With this approach, the expec- 
tation number of false detections is 1 for the entire survey. 
Applying the B07 method we find the following PFDs: 

alently pixel values of interferometric maps should follow Gaussian 
distribution. 



RT 19920826 (log(PFD)=-5.02); RT 19970205 (-2.74); 
RT 19970528 (-2.77); RT 19990504 (-4.61). With this 
more refined approach only RT 19920826 and 1 9990504 
survive. However, for reasons discussed in £|3.7.3l we have 
misgivings about RT 1990504. 

An entirely different approach^ (and in some ways or- 
thogonal to the above SNR based approach) is to look 
at the angular distribution of the transient sources with 
respect to the primary axifQ. Basic interferometry the- 
ory informs us that the dirty image is simply the Fourier 
transform of the visibility data. As such the radiomet- 
ric noise in the dirty image should be independent of 
the angular offset from the phase center. In contrast, 
the point source sensitivity decreases as one goes away 
from the pointing center and this is governed by the pri- 
mary beam response (assuming that the spectral resolu- 
tion of the survey is high enough that the delay beam 
is larger than the primary beam). Thus, once the min- 
imum SNR for detection is fixed, cosmic sources should 
be concentrated towards the pointing direction whereas 
noise spikes (masquerading as threshold point sources) 
should be uniformly distributed. 

In $B] we derive the expected distribution of cos- 
mic sources as a function of the angular offset. In 
Figure [5] we plot the expected cumulative distribu- 
tion and also the angular offset of the four sources 
which are not artifacts but whose SNR seems to be un- 
d er disp ute, namely R T 1992 0826 (EZH), RT 19 97020 5 
GEZII), RT 19970528 ( 357731) and RT 19990504 ( §57775)) ■ 
From this Figure one can see that only RT 19920826 lies 
in the expected region whereas the remaining three are 
collectively improbable. 




Fig. 5. — The cumulative probability of finding a cosmic 
source (but integrated from angular offset of infinity to zero) 
as a function of the angular offset with respect to the point- 
ing center. The angular offset is normalized in units of 
#fwhm- The sources are represented by dots and are (from 
left to right): RT 19920826, RT 19970528, RT 19970205, and 
RT 19990504. The points are deliberately placed at cumu- 
lative probability of 1%. RT 19920826 is firmly within the 
region where one naturally expects cosmic sources. 

In summary, two different statistical tests, one based 
on SNR and the other making use of the spatial signa- 

9 This test was recommended to us by J. Condon. 

10 We assume that all antennas are pointed in the same direction 
and this direction is both the pointing axis as well as the phase 
center. 
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ture provided by the primary beam, suggest that of the 
remaining four sources detected at threshold, only one, 
namely, RT 19920826 is a good detection. Thus a sim- 
ple interpretation of our re- analysis is that the rate of 
B07 transients is considerably lower than that reported 
by B07, perhaps an order of magnitude smaller. 

We acknowledge that the discrepancy between the 
analysis presented in B07 and the analysis presented here 
(and re-in vestigated) is disturbing. In the previous sec- 
tion ( §3.7p we investigate the reasons for the discrepancy 
in the SNRs and variously find possible and plausible 
causes: nagging of data, the choice of data reduction 
package (AIPS versus Miriad) and the specific method 
used to compute SNRs. However, none of these expla- 
nations are satisfactory. We are continuing this investi- 
gation but at the present time we consider this topic to 
beyond the scope of the paper. 

In contrast, the sources which we find as artifacts have 
ready explanations (see Table [2]). One source is a result 
of the file header containing a pointing direction of the 
previous pointing. Another is due to a systematic associ- 
ated with local signals (RFI). These signals do not have 
the natural fringe rate of cosmic sources and appear as 
candidates close to the phase center. Two are side-lobes 
of a stronger sources E3 Suffice to say that such rare er- 
rors will be found if one inspects sufficiently large number 
of beams! 

5. A REVISED LOOK AT THE TRANSIENT RADIO 

SKY 

In Figure |6] we plot the areal densities of three 
known transients (SN1998bw-like, type-II RSN and Swift 
J1644+57-like) and that of two expected classes (NS-NS 
mergers and orphan afterglows of long duration GRBs) 
of transients. The areal density of these five classes is 
also summarized in Table [TJ We briefly discuss each of 
these five classes of transients below. 

The traditional extra-galactic radio transients are type 
II radio supernovae. The rate that we present here is 
based on the single radio SN d etected in a blind radio 
survey by lGal-YameTafl (l2006h an d it agrees with an in- 
dependent estimate bv lLien et al.l (120111) . As illust rated 
by the example presented in iGal-Yam et all (|2006l ) . the 
main advantage in the search for radio SNe in a blind 
survey is the unique view that it provides to the oth- 
erwise hidden population of heavily obscured SNe. The 
observed rate of other types of radio SNe (e.g. ordinary 
type Ib/c SNe) are considerably lower than that of type 
II supernovae and are not discussed further. 

Radio emission is expected from both classical long 
duration gamma-ray bursts (GRBs) as well as the 
more abundant but less luminous GRBs exemplified by 
GRB 980425 associa ted with the energetic SN Ic su- 
pernova, SN 1998bw (|Galama et al.l lT998: Kul karni et al.l 
1998). The radio emission is far brighter than that of or- 
dinary core collapse SN (II, lb, Ic) and the increased vol- 
ume makes up for the intrinsically smaller birth rate. The 
recent discovery of th e energetic supernova SN2009bb 
(jSoderberg et al.ll2010| ) demonstrates that radio surveys 

11 Unfortunately, side-lobes are the exception to the expectation 
of Gaussian statistics for interferometric images. It is said that 
"the Central Limit theorem covers a large number of sins but not 
all sins." 



can find such sources without resorting to high energy 
( gamma-ray) mis s ions. 



Levinson et al.l (|2002l ) estimated the number of after- 



glows from classical GRBs and whose explosion axis is 
directed away from us ("orphan" afterglows). The ex- 
pected rate depends strongly on the poorly constrained 
7-ray beaming. On one hand this makes any rate predic- 
tion uncertain. On the other hand, even a non-detection 
by the kind of survey that we discuss below will provide 
an independent and unique constraint on the average 
opening angle of long ga mma-ray bursts, their true rate 
and total energy output dRossi et ail 120081 : iNakar et al.l 
I2002t iTotani fc Panaitescull2002D . The areal density in 
Figure [6] is derived using a typical opening angle of 10 
degrees. 

A surprising and apparently an important develop- 
ment in the field of radio transients took place just this 
year with the discovery of a radio transient associated 
with the nucleus of a modest size galaxy. The source, 
Swift J164 4+57 was initially d etected as a hard X-ray 
transient ()Burrows et al.l [20TTI ). Subsequent follow up 
found a bright, compact self-absorbed radio counterpart, 
localized at the cente r of a normal galax y at z = 0.354 
(jZauderer et al.ll201lt ITevan et al.ll20Tl . The current 
view is that Swift J1644+57 arose from a relativistic jet 
produced when a star was tidally disrupted as it passed 
too close to an oth erwise dormant su per-massive nu- 
clear black hole (e.g. iBloom et al1l2011| ). Shortly, there- 
after, a second candidate non-ther mal tidal disruptio n 
event (TDE) was recently proposed (jCenko et al.ll201ll ). 
Events such as these give us an opportunity to study 
the activity of 10 7 -10 8 M Q supermassive black holes in 
otherwise normal galaxies. 

The areal density in Figure [6] is calculated assuming 
an observed rate of 0.2 yr _1 Swift J1644 +57-like events 
and a gamma-ray beamin g factor of 10 3 ( Zau derer et al.l 
l20Tlt IBloom et al.l |20TTI) . Nominally, Swift J1644+57- 
like sources appear to be the most frequent extra-galactic 
transients that will be found in radio transient searches. 
We acknowledge that the uncertainty of both the ob- 
served rate and the gamma-ray beaming is high and the 
true rate may be significantly different 

Now we come to the most uncertain as well as poten- 
tially the most important extra-galactic radio transient 
- the merger of two neutron stars (or a black hole and 
a neutron star). It is generally accepted (or expected) 
that sho rt hard bursts are on-axis explosions of these 
mergers (jNakarl 120071 : iMetzger fc Bergerl 1201 11 1. As in 
long duration GRBs, radio emission is expected by after- 
glow (on-axis or orphan) . The rates are highly uncertain 
because there are very few observations of short hard 
GRBs. Thus there still continues to be a debate about 
the geometry of these explosions ("jetted" or not). Next, 
while the expected radio emission is straightforward to 
estimate (subject to the usual parametric uncertainties 
of the energy fractions of relativistic electrons and mag- 
netic field) an additional uncertainty is the density of the 
ambient gas (which is necessary for the production of the 
afterglow emission) . 

12 Estimates bas ed on theoretically predic te d TDE rates 
and luminosities {G iannios & Mctzgcr 2011; Bower 120111 ; 
van Velzen et aT] 120111) result in areal densities that vary by 
three orders of magnitude. The rate that we predict here is 
consistent with the upper range of these predictions. 
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Fig. 6. — Cumulative areal density of transients as a function of peak flux density for all major transient surveys. Most of the 
surveys are upper limits and the allowed phase space is above and to the right of the L-shaped symbol. The three dark blue L's 
(annotated as B2 for the 2-m onth transients from B07, Bl for the single-epoch transients from B07 and O for the lone transient 
reported in lOfek et aT]|20Tll ') are the upper limits derived as a result of the analysis presented here. These limits were derived 
b y assuming no dete ction (whence a Poisson upper limit of 3 at the 95% confidence level; see ijCj) and survey areas summarized 
in lOfek et all pHOh . 



Regardless of the uncertainty whether neutron star 
mergers are the sources of short GRBs or not, a substan- 
tial sub- and mildly relativis tic outflow is ex p ected to be 
ejected during the merger. iNakar fc Piranl ([201 ll ) esti- 
mate radio emission from these outflows. The areal den- 
sity in Figure [6] is calculated based on their estimates^, 
assuming a NS-NS merger rate of 300 Gpc _3 yr _1 and 
that any merger ejects 10 50 erg of a mildly relativistic 
outflow. We note that INakar fc Piranl (|2011l ) suggested 
that RT 19870422 was the radio emission from the re- 
mains of a neutron star merger. However, as noted in ^3] 
this source is an artifact. 

6. WAY FORWARD: NEW SURVEYS 

There are sound reasons to continue the exploration 
of the dynamic radio sky. Radio searches are an ideal 
way to discover core-collapse supernovae embedded in or 
behind dusty regions. The discovery of SN 2009bb shows 
that large radio searches can find urgently needed addi- 
tional examples of nearby low luminosity GRBs. Next, 
the many rewards of radio follow up observations of 

13 The rate density of such mergers is poorly constr ained. It 
ranges between tO to 1 4 Gp c ~ 3 yr" 1 for NS- NS mergers (Phinncy 
19911: INaravan et~all 119911 : IKalogera et~aTl 12001 lAbadie et aLI 
20171 



Swift 1644+57 (accurate localization, energetics, beam- 
ing, outflow velocity) show the tremendous diagnostic 
power of radio observations of this entirely new class of 
extragalactic transients. 

As exciting as these developments are, the search 
for new classes of radio transients has involved sev- 
eral false starts. The euphoria that followed the 
discovery of a highly dispersed (and therefore ar- 
gued to be of extr agalactic origin) millisecond burst 
(jLorimer et al.|[2007l ) was rapidly diminished by the dis- 
covery of many such bursts, pre sumably of terrestrial 
origin (|Burke-Spolaor et al.ll201ll ); but see iKeane et al.l 
(120111). Simila r ly, a long- duration transie n t foun d by 
ILevinson et~aT1 ([2002D and IGal-Yam et~aT1 (|2006D was 
later tr aced to a glitch i n the VLA on-line data taking 
system (jOfek et al.ll2010f ). ap parently affecting 0.29% o f 
all FIRST survey pointings ( Thvagaraian et aT] 120111) . 
Finally, our re-analysis (see shows that the claim 

of late time rad i o emi ssion from neutron star coalescence 
INakar fc Piranl ([201 11) is premature. 

Here we focused on the potentially new class of long 
duration radio transients reported in B07. We rule out 
six of the ten transients and cast doubts on some of the 
remaining ones. However, even with one or two survivors 
the long duration transients of B07 remain of great inter- 
est. First, even with a diminished number of transients 
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Fig. 7. — Phase space diagram showing the predicted radio lu- 
minosity versus evolutionary time scales for several types of long 
duration radio transient populations. Transparent zones indicate 
source populations which are typically optically thin, while grey 
zones indicate source populations that are expected to be optically 
thick before maximum light, evolving to an optically thin phase at 
later times . A si milar optical version of this figure can be found in 
IHau et all H2009I) and a more comprehensive figure which includes 
short duration events is in SKA Memo 97. 



the implied areal rate of B07 transients would be compa- 
rable to the recently established class of Swift J1644+57 
transients. However, unlike any other long duration ra- 
dio transient (to wit, supernovae; active stars; tidal dis- 
ruption events; gamma-ray burst afterglows, beamed or 
otherwise) the B07 transients are remarkable for the ab- 
sence of a quiescent optical counterpart. 

The event RT 19920826 survives two independent tests. 
As such it useful to speculate on the origin of this tran- 
sient. The absence of an optical (B07) and near IR qui- 
escent source could mean one of two origins. The event 
is extra-galactic in origin but the host galaxy is faint 
enough not to have been detected (as does happen for a 
few percent of long duration GRB host galaxies) or that 
it is offset from a host galaxy (as in the case for a few 
short hard bursts; [Bergcr 2009). Alternatively, the event 
is a Galactic neutron star and we have to then appeal to 
an optically invis ible Galactic (neutron star) population 
(|Ofek et al.M20T0h . 

Clearly, a new survey which can net a dozen of such 
sources (but brighter!) would help resolve the origin. For 
the Galactic hypothesis one expect no quiescent counter- 
part even at HST sensitivity whereas detectable galaxies, 
in the majority, are expected for the extra-galactic hy- 
pothesis. 

There are great gains in the discovery of new classes 
of radio transients but at the same time the path to 
true success is littered with false starts or easy specu- 
lations. The way forward must incorporate the lessons 
learnt from the false starts. We elaborate on this conclu- 
sion below. 

To start with we believe that the field of radio tran- 
sients (at least in the decimeter band) is sufficiently ma- 
ture that any new survey which just sets an upper limit 
relative to the known population is of marginal value. 
Future surveys have to be sufficiently deep and cover 
large enough sky that success (i.e. detection of a few to 
many transients) is assured. In our opinion, this means 



that a survey should be designed to find at least one or 
more Swift J1644+57-like transient (Figure [6]). 

Next, timely and multi- wavelength follow-up is es- 
sential. For examp le , the transients r eport ed in 
[ Gal-Yam etall (l200l. [Ch^gorv fc Tavlorl (fl98l and 
iBannister et al.l (|2011aM bl) have plausible origins as su- 
pernovae and Swift J1644+57-like sources (see Figure 15]). 
However, the lack of timely follow-up or deep multi- 
wavelength followup of these events preclude us from 
coming to a definitive conclusion. 

Finally, the search should be restricted to sources with 
a high level of significance. This certainly means paying 
attention to the large number of beams searched. How- 
ever, at low thresholds and with large number of beams 
(cf fj3| it would be prudent to set thresholds beyond mere 
statistical consideration^. A threshold of 9 or even 10cr 
may be appropriate. Alternatively, an immediate verifi- 
cation of a transient by deeper observation or a confir- 
mation by observations at other wavelengths would allow 
detection of transients closer to threshold. 

We start with a discussion of two recent developments. 
iBell et all (|201lD undertook an ambitious program sim- 
ilar in spirit to B07, namely the investigation of fields 
surrounding VLA calibrators. Sources with duration be- 
tween 4 and 40 d were searched for. The total integration 
time was 435 hr. No transient source in the GHz range 
with flux greater than 8 mJy was found. The authors 
place an upper limit to the areal density of 0.032 deg -2 . 
Assuming S~ 3 / 2 scaling this areal density is 4.4 deg -2 
and is well above the B07 rate. 

The FIRST survey imaged the sky in an hexagonal 
grid, in which each position in t he survey footprin t 
was observed on av e rage 3-4 times (jBecker et al.lll995T ). 
iThvagarajan et all ([201 If ) used this fact to construct 
light curves of sources detected in individual FIRST sur- 
vey snapshots. They identified 1627 variable candidates 
with variability exceeding 5er. This effort is probably the 
largest variable and transient survey ever carried out. 
One disadvantage of such a survey for transient identifi- 
cation is that the co-added images are not much deeper 
than a single epoch image. This make it hard to tell if 
an apparent transient source is really a transient or just 
a variable source that exceeded the detection threshold 
in one of the epochs. 

The limitations discussed above lead us to suggest 
a new EVLA survey specifically tailored to systemati- 
cally explore the radio sky. Following that we review 
a far more ambitious survey - the VAST key project on 
ASKAP. For the discussion below we will adopt the rates 
summarized in Table [1] The rates are specified to a flux 
density of 0.3 mJy and are extrapolated to higher flux 
densities as N(> S) oc S"" 3/2 . 

6.1. EVLA Survey 

A moderately ambitious survey with the EVLA can re- 
sult in great progress. This survey has two virtues. One, 
the EVLA offers excellent spatial resolution. Next, the 
EVLA is fully commissioned and is working to specifica- 
tions. 

14 Separately we caution that the discussion of statistics assumes 
that the underlying statistics are Gaussian to a very high degree 
of precision. As noted in the text, the sidelobes of strong sources 
add an additional source of non-Gaussian noise. 
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Specifically consider a 100 square degree survey under- 
taken in the 2-4 GHz band. An integration time of 85 s 
results in a sensitivity of 0.3 mJy (10-er). A single epoch 
covering 100 square degrees would require 50 hours. As 
can be seen from Figure [7] (taken from Table [IJ such a 
survey would have to explore a variety of time scales to 
probe the emerging classes of transient. Fifteen epochs 
could reasonably cover the range of a week to years. 

Noting the great importance of multi- wavelength imag- 
ing data, such a survey would sensibly focus on regions of 
sky where considerable multi- wavelength data (including 
radio) exists. One such region is, for example, the SDSS 
equatorial stripe (jHodge et al.ll20l"TI ). Furthermore, the 
high instantaneous sensitivity of the EVLA makes rapid 
follow up of newly transients possible. 

After the first four epochs the reference images will be 
twice as deep as the survey field. A single new epoch 
would then yield about ten Swift J1644+57-like sources 
and four supernovae. Ten such images may find a new 
example of an SN 1998bw-like event, several clear exam- 
ples of orphan afterglows and have an excellent chance 
of finding the first examples of neutron star mergers. We 
note that these different classes of objects have different 
characteristics, both in duration (Figure [7} and also in 
host magnitudes and location with respect to host galaxy. 
Therefore, it is possible to distinguish between different 
types of objects. 

As note earlier, rapid verification of a transient (either 
by additional and deeper radio observations or observa- 
tions at other wavelengths) can reduce the requirement 
for a high detection threshold. This would then require 
close rapid reduction - well within the reach of modern 
computers. 

6.2. VAST (ASKAP) 

The Variable & Slow Transient (VAST) is an approved 
key project of ASKAF0 The VAST- Wide survey aims 
to survey in the 1.2 GHz band about 10,000 square de- 
grees every day for 2 years. With 40-s integrations the 



expected single-epoch rms is 0.5 mJy (VAST Memo#l). 
Since the survey is planned to be undertaken daily the 
reference image will be built up quite rapidly. 

A comparison of data obtained at a new epoch when 
compared with the built-up reference image can be ex- 
pected to detect about six Swift J1644+57-like sources. 
Type II SN and Swift J1644+57-like sources are relatively 
long lived and so one could consider averaging the daily 
images (to say 10 days). The resulting summed data has 
a sensitivity of 0.16 mJy. The expected number of tran- 
sients per such summed image is nominally 32 (Type II 
SN), 82 (Swift J1644+57), 8 (orphan afterglows) and 4 
(neutron star coalescences). 

To sum up, the dynamic radio sky remains a rich area 
for exploration. Based on what we know about the areal 
rates for the known transient sources, future synoptic 
radio imaging surveys are expected to yield substantial 
numbers of exotic transients. Such surveys will also pro- 
vide the definitive test for the B07 population. 
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APPENDIX 

A. PROBABILITY DENSITY FUNCTION OF MTH MAXIMUM 

In this section, our goal is to compute the probability density function of the mth highest value of hj, j = 1, 2, 3...n. 

Let p(h) be the probability density function of hj with P(h) = ^ oo p{h)dh being the cumulative function. We denote 

the m th highest value by H m . Thus the maximum of the series of measurements is H\ and the minimum value is H n . 

Let p{H m ) be the probability density function of h — H m . This means that at least one of the measurements lies is 
in the range [H m , H rn + dH rn ]. The probability density for this event is p(H rn ). Next, then n — m measurements must 
lie below this range and m — 1 above this range. The probability for any value to be smaller than H m is P(H m ) and 
the probability for a value to higher than H m is is 1 — P(H m ). This is now a binomial distribution with n — 1 total 
values. Thus the probability density function for H rn is 

m — 1 

(Al) 

The first combinatorial factor of n accounts for the possibility that H m can occupy any position in the sequence. The 
second combinatorial factor, " _1 C m _i account for the combinations satisfying the condition that n — m values lie 
below H m and m — 1 lie above H m . For both the maximum (to = 1) and minimum (m = n) Equation IA1I simplifies 
to that expected from basic considerations. 
Now let us consider the specific case where h follows Gaussian statistics: 

P (h) = -}=eM~h 2 /2) (A2) 

V 27T 

where h is normalized in units of a. The probability that an event is extreme or lies within the range ±/i is 

p(h)dh = ^ exp(-z 2 )dz = erf (h/V2). (A3) 

-h V 71 " Jo 

The probability that an event is extreme in only one sense, maximum or minimum, and lies outside the range [— oo, h] 
(say) is thus 

PW = i[l + erf(A) 

where erfc(a;) = 1 — erf(x). 



p{H m ) = np(H m ) x 



in, — mlllm, — 1 ) 



■P(H m ) n ~ m l-P(H m ) 



2 erfc (^ 



(A4) 
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Consider the case where n» 1 (say 10 6 or more) and m is small, say 10. Then we can approximate n — 1 ~ n, 
n — 2 w n, n — to + 1 ss n. Furthermore using the approximation, (1 — x/n) n exp(— x), we find in the limiting 
case where H m is greater than a few [so that 1 — P{H m ) <C 1]: 



p(H„ 



(to- 1)! 



p(H m ) exp 



(n 



5-0) 



(A5) 



B. DISTRIBUTION OF SOURCES WITHIN THE PRIMARY BEAM 



Provided that there is sufficient spectral resolution, the response of an interferometer to a source follows the antenna 
response function ("primary beam")0 We will assume that the response function is azimuthally symmetric and 
specified by g(9) where 9 is the angular offset from the pointing axis. 

Let the areal density of sources with flux density greater than S be a power law, say, N(> S) oc S a . For Euclidean 
geometry and most reasonable luminosity functions, a — —3/2. Next, we note that detection is really finding sources 
at a given SNR and above. Fortunately the noise distribution for an interferometric image is uniform. Thus a source 
with a given flux density will have an SNR, S, that scales with the primary beam response, S oc Sg(9). The number 
of sources above a certain SNR and contained outside an angular radius of 9q is 



n(>S;>0o)« / ^ei-^-YdO 



For the specific case of a Gaussian beam, 



(Bl) 



(B2) 



where the tra ditio nal "full width at half maximum" (FWHM) is 6*fwhm 
into Equation IB II we obtain 

n{>S->6 Q )K (s/g(9 )Y 



^/ln(256)#*. Substituting Equation IB2I 

(B3) 



Half the sources will be detected outside the radius 9h = y / — ln(4)/a 9*. For a = —3/2 we obtain 9h = -\/ln(16)/3 0* ~ 
0.970* . The expression 9h = yjl /6 $fwhm ~ 0.4 #fwhm is more useful. A plot of n(> S; > 9) can be found in Figure[8l 




Fig. 8. — The expected distribution of sources above a threshold SNR as a function of 6/8fwhm from the pointing aixs and 
assuming a — —3/2. The response curve is normalized by insisting that the integral of the curve (from 6 = to 9 = oo) is unity. 
Fifty percent of the sources are within 0.416*pwhm (marked by a pentagram) and 97% within 0.92#fwhm (marked by a cross; 
at this radius the primary beam gain is only 0.1 relative to that on-axis). 



C. UPPER LIMIT FOR A NON DETECTION (POISSON) 

It is not unusual to find no source after undertaking a survey. We wish to determine an upper limit to the number 
of sources that were being searched. The number detected is given by Poisson statistics. The probability of finding r 
sources is then given by 

p(r) = ^ exp(-A) (CI) 

16 We assume that the phase center coincides with the pointing 
axis of the primary antenna. 
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where A is the Poisson parameter and equal to (r) . 

Our goal is to determine the maximum value of A given a non-detection. As the value of A is increased the probability 
of detection, by which we mean the probability of detecting one or more events, also increases. This probability is 
p(l) + p(2) + ... which we note is 1 — p(0). This probability can be set to the desired confidence level, V and thence 

V = 1 -p(0) = 1 - exp(-A). (C2) 

The reader with a stronger physical bent may find the complement, p(0) =1 — V, more appealing. Regardless, we find 
A = [3, 4.6, 6.9] at a confidence level of [95%, 99%, 99.9%]. 

D. NUMBER OF INDEPENDENT BEAMS 

Here we compu te n, the number of independent beams for the VLA data that went into the analysis of B07 and 
IQfek et all (l20ll . For B07, a circular region with a radius of two times the half-power radius was searched for each 
epoch. The number of independent beams per epoch is the ratio of that area to the area of the synthesized beam in that 
particular epoch. For individual epochs, this value ranged from as small as 10 3 to as large as 10 6 . The total numbers 
of independent beams for the 5 and 8.4 GHz data are 9.3 x 10 7 and 4.5 x 10 7 , respectively. The smaller number of 
independent beams for 8.4 GHz are due to the tapering of the visibility data, increasing the typical synthesized beam 
size a t 8.4 GHz. 

For lOfek etahl (|2011l ) the search was made for transients in 4.65'radius circular region. The FWHM of the synthesized 
beam is « 4 arc-seconds. Though the goal was 16 epochs per pointing we only achieved an average of 15.7 images. 
With 141 epochs we derive n = 1.1 x 10 7 . 



